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GROUND- WAVE  PATH  LOSS  OVER 
ICE-COVERED  SEA  AT  FREQUENCIES 
BETWEEN  0.1  AND  10  MHz 


I.  INTRODUCTION 

Electromagnetic  energy  may  travel  between  two  terminals  by  any  of  several  propagation 
mechanisms  depending  on  the  distance  and  nature  of  the  path  separating  the  terminals.  Some  propa¬ 
gation  modes  such  as  meteor  scatter  are  found  to  vary  on  a  short-term  basis  while  others,  including 
sky-wave  ionospheric  refraction,  may  vary  in  a  longer,  yet  still  time-dependent  manner. 

When  the  distance  between  terminals  is  within  or  somewhat  beyond  line-of-sight  (LOS),  the 
ground-wave  mode  may  be  the  predominant  propagation  mechanism.  In  general,  the  ground  wave 
consists  of  direct  and  ground-reflected  components  (existing  when  the  terminals  are  within  LOS),  and 
a  surface-wave  component  that  may  persist  to  distances  well  beyond  LOS  depending  on  the  frequency 
and  ground  characteristics. 

The  nature  of  the  ground  wave  has  been  studied  and  reported  by  numerous  investigators  for  a 
variety  of  path  characteristics  including  flat  and  spherical  smooth  earth.  While  a  number  of  useful 
paths  have  been  studied,  including  seawater  and  lossy  earth,  in  most  cases  the  electrical  parameters  of 
the  ground  have  been  assumed  to  be  uniform  with  depth.  Most  paths  of  practical  interest  including 
the  important  case  of  intratask  force  communication  in  the  open  sea  can  be  analyzed  using  these 
models.  Barrick  [1]  provides  a  useful  extension  of  the  theory  for  a  rough  sea. 

The  variation  of  the  ground  wave  excited  over  a  flat  earth  surface  whose  electrical  parameters 
are  not  uniform  with  depth  was  studied  by  Wait  [2]  and  extended  to  the  spherical  surface  case  [3], 
With  this  analysis,  it  became  possible  to  predict  the  important  path  characteristics  for  many  links  of 
current  interest  to  the  Navy,  including  propagation  over  an  ice-covered  sea.  In  papers  by  Hill  and 
Wait  [4,5],  it  was  shown  that  at  high  frequency  (HF),  the  effect  of  the  ice  layer  is  to  produce  an 
enhanced  field  strength  (over  the  seawater  only  case)  at  shorter  ranges,  but  with  serious  degradation 
as  the  distance  is  increased  beyond  a  point  that  is  strongly  dependent  on  frequency. 

Whether  propagation  occurs  over  homogeneous  or  stratified  earth,  the  ground-wave  mode 
exhibits  a  number  of  important  advantages  over  other  propagation  types.  These  include: 

•  path  loss  is  not  time  dependent  over  intervals  associated  with  typical  message  traffic; 

•  optimum  frequency  can  be  predicted  in  advance  of  communication  sessions  and  does  not 
require  active  channel  probing  as  may  be  the  case  for  sky-wave  or  meteor-scatter  modes; 

•  frequency  dispersion  is  insignificant  (although  for  propagation  over  sea  ice,  this  may  be  a  fac¬ 
tor  at  some  ranges  and  frequencies); 

•  unaffected  by  natural  or  man-made  ionospheric  disturbances,  and 
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•  there  is  a  built-in  advantage  over  interceptors  and  jammers  located  beyond  the  communication 
range  because  of  the  sharp  increase  in  attenuation  at  such  distances. 

The  purpose  of  this  report  is  to  extend  the  work  presented  by  previous  investigators  to  include 
the  low-frequency  (LF)  and  medium-frequency  (MF)  range  beginning  at  0.1  MHz  where  the  electrical 
properties  of  an  ice-covered  sea  will  be  shown  to  produce  enhanced  link  performance  over  distances 
of  practical  importance. 

II.  GROUND- WAVE  ATTENUATION  FUNCTION 

Under  unrestricted  conditions  and  following  the  general  analysis  of  Hill  and  Wait  [6],  the  elec¬ 
tric  field  strength,  E ,  produced  at  a  distance  from  a  vertical  electric  dipole,  is  given  by: 

E=E0W ,  (1) 

where  E0  is  the  field  strength  produced  by  the  same  source  over  a  plane,  perfectly  conducting  surface 
and  W  is  the  attenuation  function  that  includes  effects  caused  by  the  spheroidal  and  stratified  surface 
shown  in  Fig.  1.  Medium  1  is  the  layer  having  a  conductivity  ox  and  dielectric  constant  ej.  The 
lower  region,  assumed  to  extend  to  great  depth,  is  characterized  by  conductivity  o2  and  dielectric  con¬ 
stant  62*  The  upper  half-space  is  assumed  to  be  equivalent  to  free  space  with  no  conductivity  and 
dielectric  constant  e0.  The  permeability  in  all  regions  is  equal  to  fi0. 


Vertical  Electric 


Figure  1  —  Propagation  over  stratified,  spherical  earth 


Electric  Field 
Strength 


When  the  height  of  the  dipole  source  and  the  point  at  which  the  field  strength  is  measured  are 
both  at  the  upper  surface  of  medium  1,  W  can  be  written  as  the  residue  series: 
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The  surface  impedance  Z  at  the  boundary  between  medium  1  and  the  upper  half-space  is  given 
[41  by: 
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The  values  ts,  which  make  up  the  poles  of  Eq.  (2)  and  appear  in  the  definition  for  X,,  are  roots 
of  the  equation. 


where 


dw(t) 

dt 


-  qw(t)  =  0, 


(4) 


w(t)  =  Vtt  [Bi  (t)  —  iAi  (r)] . 


The  symbols  Ai{t)  and  Bi(t)  refer  to  Airy  functions  [7],  The  system  of  Eqs.  (3)  and  (4)  can  be 
solved  iteratively  for  the  roots  ts  (the  program  listing  given  in  Appendix  A  provides  this  procedure). 
However,  as  pointed  out  by  Hill  and  Wait  [4],  for  the  case  of  Arctic  ice  overlaying  a  highly  conduc¬ 
tive  sea,  a  reasonable  approximation  for  Z  in  Eq.  (3)  can  be  obtained  by  setting  \s  =  0. 

A  general  procedure  for  solving  Eq.  (4)  for  the  roots  ts  is  given  in  Hill  and  Wait  [6]  and  is 
briefly  summarized  in  Appendix  B.  Two  simpler  forms  for  the  ground-wave  attenuation  function  W 
can  be  used  at  relatively  short  numerical  distances  x .  The  first  of  these,  called  the  small  curvature 
expansion,  is  a  modification  to  the  flat  earth  attenuation  function  for  a  spherical  surface  and  is  useful 
for  large  |  q  | .  The  second  form,  called  the  power  series  representation,  is  most  useful  at  short 
ranges  and  small  |  q  |  (e.g.,  a  nonlayered,  highly  conducting  surface  such  as  open  ocean).  Both  of 
these  approximations  along  with  their  regions  of  validity  are  presented  in  Appendix  B. 

III.  NUMERICAL  RESULTS 

The  equations  for  the  ground-wave  attenuation  function  have  been  solved  for  the  case  of  sea  ice 
of  various  thicknesses  overlaying  seawater.  In  the  results  that  follow,  a  dielectric  constant  er2  =  80 
and  conductivity  o2  =  4  Siemans  was  assumed  for  the  seawater  constituting  medium  2.  (A  descrip¬ 
tion  of  a  special,  multilayer  case  follows  in  which  a  lower  conductivity  seawater  layer  exists  just 
below  the  ice  sheet  as  a  result  of  melting.)  For  the  sea-ice  layer,  the  conductivity  and  dielectric  con¬ 
stant  are  known  to  be  dependent  on  frequency,  temperature,  and  the  age  of  the  ice  with  cold  or  mul¬ 
tiyear  ice  being  much  less  conductive  than  warm  or  first-year  ice.  Since  the  variation  of  conductivity 
with  frequency  is  less  pronounced,  a  constant  value  for  ol  =  0.000333  Siemans  was  used  for  the 
sea-ice  layer  of  medium  1.  This  corresponds  to  one-year-old  ice  at  -10°C  as  given  by  Wentworth 
and  Cohn  [8],  The  variation  of  dielectric  constant  er]  with  frequency  is  more  strongly  dependent  on 
frequency  and  is  also  derived  from  data  given  in  Ref.  8  (see  Fig.  2).  Another  useful  quantity  that  can 
be  used  to  describe  the  electrical  quality  of  the  ice  sheet  is  the  loss  tangent  defined  as  tan  5  =  o/a>ere0 


One  of  the  major  difficulties  in  analyzing  a  sea-ice  path  is  the  considerable  range  of  electrical 
parameters  that  is  possible  for  the  ice  layer  depending  on  season  and  geographic  location.  It  is  impor¬ 
tant  to  note  that  the  results  presented  here  are  based  on  the  particular  values  chosen  for  the  dielectric 
constant  and  conductivity  of  the  ice  layer.  Although  they  are  believed  to  be  representative  of  a  typi¬ 
cal  Arctic  path,  some  variation  in  link  characteristics  can  be  expected  depending  on  the  specific  geog¬ 
raphy  chosen. 

The  total  path  loss  is  presented  in  the  results  that  follow,  including  effects  due  to  the  attenuation 
function  W  as  well  as  the  distance  dependence  of  the  term  E0  in  Eq.  (1).  Furthermore,  it  is  assumed 
that  for  any  given  case,  the  electrical  properties  of  the  surface  over  the  complete  path  are  constant  and 
that  the  thickness  of  the  layer  is  uniform.  Although  techniques  are  available  for  analysis  of  a  mixed 
path  in  which  the  ice  layer  varies  in  its  electrical  properties  and  thickness,  this  has  not  been  con¬ 
sidered  here. 
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Figure  2  —  Relative  dielectric  constant  variation  as  a  function 
of  frequency  used  in  calculations  of  ground-wave  path  loss 


Figures  3  to  12  show  the  numerical  results  of  the  path  attenuation  in  decibels  as  a  function  of 
distance  in  kilometers  for  frequencies  between  0.1  and  10  MHz.  Table  1  provides  a  summary  of  the 
data  presented  in  each  figure  as  well  as  the  dielectric  constant  and  loss  tangent  used  in  the  calculation. 
Each  figure  also  includes  two  curves  for  a  uniform  earth  (medium  1  =  medium  2),  one  for  seawater 
only  and  one  for  sea  ice  only,  for  comparison  with  the  layered  cases.  Ice  thicknesses  up  to  7  m  were 
chosen  for  study  since  thicker  layers  are  rarely  encountered  over  a  large  horizontal  scale.  A  typical 
thickness  of  3  m  is  frequently  given  for  the  central  Arctic.  Only  thinner  ice  sheets  are  considered  at 
the  higher  frequencies  because  of  the  much  stronger  effect  on  the  path  attenuation. 


Figure  3  —  Ground-wave  attenuation  as  a  function  of  distance  over  smooth  spherical  earth 
for  uniform  ice  and  seawater  surfaces  and  for  three  layered  cases.  The  frequency  is  100 
kHz. 
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Figure  4  —  Ground -wave  attenuation  as  a  function  of  distance  over  smooth  spherical  earth 
for  uniform  ice  and  seawater  surfaces  and  for  three  layered  cases.  The  frequency  is  200 
kHz. 


Figure  5  —  Ground-wave  attenuation  as  a  function  of  distance  over  smooth  spherical  earth 
for  uniform  ice  and  seawater  surfaces  and  for  three  layered  cases.  The  frequency  is  300 
kHz. 
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Figure  6  —  Ground- wave  attenuation  as  a  function  of  distance  over  smooth  spherical  earth 
for  uniform  ice  and  seawater  surfaces  and  for  two  layered  cases.  The  frequency  is  500 
kHz. 


Figure  7  —  Ground-wave  attenuation  as  a  function  of  distance  over  smooth  spherical 
earth  for  uniform  ice  and  seawater  surfaces  and  for  three  layered  cases.  The 
frequency  is  700  kHz. 
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Higure  8  —  Ground-wave  attenuation  as  a  function  of  distance  over  smooth  spherical 
;arth  for  uniform  ice  and  seawater  surfaces  and  for  three  layered  cases.  The 
'requency  is  1  MHz. 


RANGE  (km) 


Figure  9  —  Ground-wave  attenuation  as  a  function  of  distance  over  smooth  spherical 
earth  for  uniform  ice  and  seawater  surfaces  and  for  two  layered  cases.  The  frequency 
is  3  MHz. 
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Figure  10  —  Ground- wave  attenuation  as  a  function  of  distance  over  smooth  spherical 
earth  for  uniform  ice  and  seawater  surfaces  and  for  two  layered  cases.  The  frequency 
is  5  MHz. 


DISTANCE  (km) 

Figure  11  Ground- wave  attenuation  as  a  function  of  distance  over  smooth  spherical 
earth  for  uniform  ice  and  seawater  surfaces  and  for  three  layered  cases.  The 
frequency  is  7  MHz. 


8 


NRL  REPORT  9041 
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Figure  12  —  Ground-wave  attenuation  as  a  function  of  distance  over  smooth  spherical 
earth  for  uniform  ice  and  seawater  surfaces  and  for  three  layered  cases.  The 
frequency  is  10  MHz. 


Table  1  —  Data  Summary 


Figure 

Frequency 

(MHz) 

Ice  Thickness 
(m) 

Relative  Dielectric 
Constant 

Tan  8 

..  3 

0.1 

3,5,7 

22  . . 

"  1.2 

4 

0.2 

3,5,7 

16 

1.9 

5 

0.3 

3,5,7 

13 

1.3 

6. 

0.5 

3,5 

10 

1.0 

7 

0.7 

3,5,7 

8 

1.0 

8 

1.0 

1,3,5 

7 

1.0 

9 

3.0 

1,3 

6 

0.333 

10 

5.0 

1,3 

6 

0.200 

11 

7.0 

.5,1,3 

6 

0.143 

12 

10.0 

.5,1,2 

6 

0.100 
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In  the  calculation  of  these  results,  one  of  the  three  techniques  described  in  Appendix  B  was 
used.  In  general,  for  a  given  case,  either  the  power  series  or  small  curvature  approximation  was  used 
at  the  short  range  and  the  residue  series  for  the  longer  ranges.  To  obtain  a  smooth  fit  at  the  transition 
range,  50  terms  were  used  in  the  residue  series  summation. 

The  general  characteristic  of  the  variation  of  path  loss  for  the  layered  case  can  be  described  in 
terms  of  three  regions.  The  lower  frequency  curves  (Figs.  3  to  6)  show  only  the  first  of  these  zones, 
a  region  at  short  ranges  where  the  surface  wave  mode  dominates  and  the  field  strength  is  enhanced 
over  the  seawater  only  case.  The  effect  is  minor  for  the  thicknesses  studied  (and  of  practical  interest) 
at  the  lowest  frequencies  but  is  significant  at  300  and  500  kHz  where  the  path  loss  is  as  much  as  25 
dB  better  than  seawater  alone  at  a  range  of  1000  km  (546  nmi).  The  region  of  enhanced  propagation 
is  evident  in  all  higher  frequency  curves  although  at  distances  that  become  shorter  as  the  frequency 
increases.  At  the  two  highest  frequencies  studied,  5  and  7  MHz,  the  enhanced  zone  still  exists  but  at 
distances  too  short  to  be  useful  for  most  communication  purposes. 

The  second,  intermediate  zone  is  a  region  in  which  the  field  strength  decreases  rapidly  and  non- 
monotonically  with  distance.  The  great  increase  in  path  loss  over  this  interval  is  caused  by  the  decay 
of  the  surface-wave  component  of  the  total  field  in  this  region.  At  some  distance  beyond  the  onset  of 
this  rapid  increase  in  path  loss,  the  surface-wave  mode  decays  to  a  magnitude  that  is  comparable  to 
the  other  normal  ground-wave  modes,  and  an  interference  region  consisting  of  maxima  and  minima  in 
the  attenuation  function  is  evident.  Although  not  shown  here,  there  is  a  corresponding  oscillatory 
effect  in  the  phase  of  the  electric  field  as  a  function  of  range  [5]  rendering  this  portion  of  the  path 
dispersive. 

Because  the  surface-wave  mode  has  a  greater  exponential  attenuation  than  the  other  modes  mak¬ 
ing  up  the  total  ground-wave  field,  its  effect  on  the  behavior  of  the  total  path  loss  is  no  longer  impor¬ 
tant  beyond  the  intermediate  region.  In  the  third  zone,  the  attenuation  function  returns  to  a  smooth, 
monotonically  decreasing  function  of  range  but  at  a  much  greater  loss  than  for  seawater  alone,  in  spite 
of  the  existence  of  a  seawater  substrate.  Moreover,  in  this  region  the  shape  of  the  layered  case 
curves  matches  the  ice-only  curve,  indicating  that  the  ground  wave  in  this  zone  is  dominated  by  the 
ice  and  not  by  the  seawater  substrate.  It  is  readily  apparent  that  at  frequencies  above  approximately  2 
to  3  MHz  for  ice  thicknesses  likely  to  be  encountered  over  deep  Arctic  ground- wave  paths,  the  path 
loss  is  much  greater  than  for  an  open  ocean  geometry,  limiting  the  range  of  available  frequencies  for 
communication  purposes  to  the  MF  and  lower  HF  ranges. 

This  effect  is  seen  more  graphically  in  Fig.  13  where  the  path  loss  is  plotted  as  a  function  of 
frequency  for  an  ice  thickness  of  3  m  for  four  distances,  50,  200,  800,  and  1600  km.  For  com¬ 
parison,  the  dashed  curve  shows  the  attenuation  as  a  function  of  frequency  for  a  smooth  seawater-only 
case  at  a  range  of  50  km.  (In  the  solid  curves,  the  oscillatory  behavior  associated  with  the  layered 
cases  at  intermediate  ranges  has  been  smoothed  for  clarity.)  Several  interesting  characteristics  are 
evident  in  Fig.  13  including: 

•  For  any  given  range,  at  the  lowest  frequencies,  the  layered  and  uniform  earth  curves  are 
asymptotic.  Apparently  the  layer  thickness  (and  its  effect)  is  trivial  at  these  frequencies  in 
comparison  with  the  wavelength  of  the  electric  field  in  air. 

•  There  is  a  frequency  or  set  of  frequencies  where  the  layered  earth  propagation  enhancement 
over  uniform  seawater  is  optimized.  For  a  distance  of  50  km,  this  frequency  is  near  1  MHz 
while  at  the  much  longer  distance,  1600  km,  the  optimum  frequency  has  dropped  to  approxi¬ 
mately  500  kHz. 

•  The  expected  increasing  trend  of  attenuation  with  frequency  that  is  evident  in  the  uniform  sea¬ 
water  and  short-range  layered  cases  is  reversed  for  the  two  longer  distance  cases  plotted.  At 
1600  km,  for  example,  frequencies  near  500  kHz  have  less  path  loss  than  lower  frequencies  in 
the  range  100  to  300  kHz. 
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Figure  13  —  Ground- wave  attenuation  over  spherical  earth  as  a  function  of  frequency 


•  Each  of  the  solid  curves  contains  a  sharp  knee  in  its  variation  with  frequency.  This  feature 
becomes  more  prominent  as  the  distance  increases.  At  1600  km,  for  example,  the  differ¬ 
ence  in  path  loss  between  600  and  800  kHz  is  on  the  order  of  40  dB.  The  effect  is  similar 
to  waveguide  cutoff  and  highlights  the  importance  of  careful  choice  of  frequency  in  a  prac¬ 
tical  case. 

•  Frequencies  in  the  HF  band  (3  to  30  MHz)  are  of  value  only  for  short  distance  links  via 
ground-wave  propagation.  Longer  distance  links  are  better  addressed  using  frequencies  in 
the  MF  range  (0.3  to  3  MHz).  The  availability  of  communication  hardware  capable  of 
matching  this  need,  however,  is  limited. 

Figure  14  shows  the  variation  of  path  loss  as  a  function  of  layer  depth  for  10  and  50  km  at  a 
fixed  frequency  of  7  MHz.  The  two  curves  approach  the  seawater-only  case  (shown  as  horizontal 
lines  at  the  ordinate  axis)  as  the  thickness  is  reduced  to  zero.  At  the  right  edge  of  the  figure,  the 
curves  asymptotically  approach  the  case  of  propagation  over  uniform  sea  ice  as  the  thickness  increases 
indefinitely.  For  small  ice  depths,  a  modest  enhancement  over  uniform  seawater  is  shown  for  both 
curves.  As  the  thickness  increases,  the  path  loss  increases  rapidly  reaching  a  maximum  attenuation  at 
a  thickness  of  approximately  4.745  m  and  then  reversing  this  trend  to  reach  a  minimum  near  9.44  m. 
A  succession  of  maxima  and  minima  follow  the  envelopes  shown  for  further  thickness  increases. 
This  effect  is  a  consequence  of  the  impedance  transformation  effect  due  to  the  ice  layer  and  described 
in  Eq.  (3).  The  underlying  highly  conductive  sea  substrate  is  alternately  transformed  by  this  mecha¬ 
nism  in  a  manner  similar  to  a  conventional  transmission  line  to  produce  a  surface  impedance  that  is, 
alternately,  a  good  conductor  and  a  poor  conductor.  Since  the  ice  is  not  a  perfect  dielectric  and 
possesses  some  loss,  successive  maxima  and  minima  do  not  reach  preceding  levels  and  the  underlying 
sea  is  eventually  masked  by  ice  losses  as  the  thickness  continues  to  increase. 

Figure  15  shows  three  cases  of  a  multilayer  geometry  at  700  kHz  for  comparison  with  a  uniform 
3-m  thick  sea-ice  layer.  In  the  curve  identified  as  case  A,  six  layers  of  sea  ice,  each  0.5-m  thick  are 
used  in  the  computer  model  to  form  a  piecewise  approximation  to  the  variation  of  conductivity  with 
depth  in  a  3-m  thick  sheet.  Table  2  gives  the  values  used  that  were  obtained  from  Wentworth  and 
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Figure  14  —  Ground-wave  attenuation  at  fixed  distances  of  10  and  50  km 
as  a  function  of  ice  layer  thickness  at  7  MHz 


Figure  15  —  Ground-wave  attenuation  as  a  function  of  distance  for  propagation  over  an 
ice  covered  sea  surface  where  the  electrical  properties  of  the  layer  vary  with  depth 
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Table  2  —  Ice  Electrical  Properties  Used  in  Fig.  15 


Depth 

Curve  A 

Curve  B 

Curve  C 

(m) 

er 

a 

a 

er 

a 

0.5 

5.5 

0.000100 

5.5 

0.000154 

5.5 

0.000100 

1.0 

6.5 

0.000159 

6.0 

0.000210 

6.5 

0.000159 

1.5 

8.0 

0.000251 

6.0 

0.000280 

8.0 

0.000251 

2.0 

10.0 

0.000398 

6.5 

0.000341 

10.0 

0.000398 

2.5 

12.0 

0.000631 

7.0 

0.000423 

12.0 

0.000631 

3.0 

18.0 

0.001000 

7.5 

0.000467 

18.0 

0.001000 

3.5 

80.0 

1.000000 

80.0 

4.000000 

80.0 

4.000000 

>4.0 

80.0 

4.000000 

80.0 

4.000000 

80.0 

4.000000 

Cohn  [10].  In  certain  seasons,  ice  melt  can  produce  a  layer  of  fresh  water  that  tends  to  stratify  at  the 
ice-water  boundary  since  it  is  less  dense  than  regular  seawater.  For  case  A,  this  fresh-water  layer  has 
been  included  in  the  model  with  a  thickness  of  0.5  m  and  a  conductivity  a  =  1  Sieman.  Curve  C  is 
identical  in  its  parameters  to  Curve  A  with  the  exception  that  the  layer  of  fresh  water  is  not  present. 
Curve  B  contains  data  obtained  from  a  core  sample  in  the  marginal  ice  zone  during  an  Arctic  field 
test. 


Differences  in  the  curves  are  not  significant  until  the  distance  increases  to  beyond  approximately 
200  km.  At  the  greatest  range  shown,  the  maximum  difference  between  curves  is  on  the  order  of  10 
dB.  Curves  A  and  C  are  plotted  for  ice  parameters  that  are  significantly  less  lossy  than  for  the  ice 
used  for  Curve  B,  but  Curve  A  shows  significantly  more  attenuation  at  the  greatest  range  because  of 
the  less  conductive  seawater  substrate.  Note  that  curves  B  and  C,  which  describe  a  more  realistic  ice 
layer,  are  still  within  3  dB  of  the  uniform  parameter  3-m  case  even  at  the  longest  range  considered. 

IV.  SUMMARY 

The  ground- wave  propagation  mode  has  been  studied  for  applications  involving  a  sea-ice  path. 
In  general,  the  advantages  of  this  propagation  mechanism  are: 

•  The  path  loss  is  not  time-dependent  over  intervals  associated  with  normal  communication  mes¬ 
sage  length.  The  optimum  frequency  can  be  predicted  in  advance  of  communication  sessions 
and  does  not  require  active  channel  probing  as  may  be  the  case  for  skywave  or  meteor  scatter 
propagation. 

•  Frequency  dispersion  is  insignificant  (although  for  propagation  over  ice,  this  may  not  be  true 
at  some  ranges  and  frequencies). 

•  There  is  a  large  margin  against  jammers  or  interceptors  who  are  located  beyond  the  communi¬ 
cation  range. 

•  Ground-wave  propagation  is  unaffected  by  natural  or  man-made  ionospheric  disturbances  and, 
in  fact,  may  be  enhanced  during  such  events  since  external  noise  and  interference  levels  can 
be  expected  to  diminish  in  response  to  ionospheric  failure. 
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In  the  Arctic,  some  of  these  advantages  are  magnified.  The  normal  knee  in  the  attenuation 
curve  for  seawater  is  changed  appreciably  for  an  ice-covered  path  as  a  result  of  the  rapid  decay  of  the 
surface  wave  component  of  the  total  field.  Correct  choice  of  operating  frequency  can  simultaneously 
optimize  communication  link  performance  while  discriminating  against  potential  interceptors  or  jam¬ 
mers. 


It  has  been  shown  that  the  path  characteristics  are  sensitive  to  the  electrical  properties  of  the  sur¬ 
face,  particularly  the  ice  thickness,  but  also  the  conductivity  and  dielectric  constant  of  the  ice.  Dur¬ 
ing  those  times  of  the  year  when  the  ice  melts,  path  loss  will  probably  be  affected  by  the  presence  of 
a  fresh  water  layer. 

The  general  behavior  of  the  attenuation  as  a  function  of  range  is  given  to  a  good  degree  of  accu¬ 
racy  by  considering  the  layer  to  be  of  uniform  electrical  characteristics  throughout  its  thickness.  An 
analysis  of  the  effect  of  changes  in  layer  thickness  over  the  path  length  has  not  been  done.  It  is 
expected  that  the  higher  frequencies  will  be  most  strongly  affected  by  ice  thickness  variability  since 
small  changes  will  represent  a  significant  portion  of  a  wavelength. 

In  the  Arctic,  the  best  frequencies  for  use  in  ground-wave  links  range  from  100  kHz  or  below  to 
approximately  2  MHz  depending  on  path  length.  When  the  attenuation  characteristics  are  combined 
with  other  system  parameters  such  as  antenna  efficiency  and  external  noise  level,  the  optimum  com¬ 
munication  frequency  is  expected  to  be  readily  identifiable  as  a  peak  in  system  performance  over  a 
small  range  of  frequencies  in  the  MF  band. 

Complete  system  analysis  must  also  consider  the  effects  of  ionospheric  refraction  on  link  perfor¬ 
mance.  Electromagnetic  energy  reaching  a  distant  terminal  by  this  mode  may  be  of  comparable  mag¬ 
nitude  to  the  ground  wave  leading  to  a  zone  of  interference  between  the  two  components.  In  addi¬ 
tion,  the  Artie  ionosphere  is  known  to  be  highly  variable  with  frequent  periods  of  ionospheric  failure. 
Under  these  conditions,  external  noise  and  interfering  signals  that  reach  the  Arctic  by  sky  wave 
modes  can  be  expected  to  decrease  markedly  leading  to  improved  link  performance. 
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Appendix  A 
PROGRAM  LISTING 


I  0 
20 
30 


OVERLAP 
RAD 

COM  U  .EG >U0  ,Er (1 0 ) ,Si (10) ,Hi< 1 0 ) ,Gr ( 1 0) ,Ci (1 0  > ,Kfs,Re,Nrt ,N1* 

i n n 


Nrt-100 
-DISP.  " 


80 

90 

100 

no 

120 


UAXX— =_  MATRIX — INITIALIZATIONS _ 


PRINTER  IS  16 
CALL  C»xp<0, -PI/3, Rr,Ri) 

1--0 1879297 , 3 . 24819758  ,  4  .  8200992.1  ,6 , 16330736  >7 . 37217726,8 . 48848673 
DATA  9 . 53544905 , 1  0 . 52766040 , 1 1  .-47505663,12.38478837 
RESTORE  100 
12LQ _ FOR  J-l  TO  in _ 

140  READ  X5  - - 

150  Yit  r  < J)«X5*Rr 

1611  Yl.l  i  CJ4«X5.*Ri _ 

170  NEXT  J  "  '  " 

Jon  08794944,3 .52055983,6.78670809,7.94413359,9. 02265085 

- DATA  10  ,-04017434.11  ■  flimiPAin  ,  1  1  P3AQ1  3SA  1  ?  .  RPR77A7"} _ 

200  RESTORE  180  - 

210  FOR  J»1  TO  10 

22H—  -READ  X5 _ _ _ 

230  Yi2r( J)-X5*Rr  ‘ 

240  Yi2i(J)«X5*Ri 

250 - NEXT  J _ _ _ 

260  FOR  J«U  TO  Nrt ”  “  "  ~  "  “  “  ~ - - - ’ - 

270  Ai*»3*PI*(4*J-l )/8 

2l  ) - Ai*p«3*PL*M4*J-30/8 _ — _ _ _ _ 

290  CALL  Ai«ro( Ais,X3>  “ 

300  CALL  Az»r op (Ai*p , Y5) 

-34-0 - Y4-1  r  ( Jl»Y5*Rr _ ’ _ 


320 
330 
-34  a 
330 
360 
370 


Yilif J)«Y5*Ri 
Yi2r( J)*X5*Rr 

Y424  CJ-)  »X3*Ri _ 

NEXT  J 
MAT. T*r-ZER 
-MAT — Ts.i*ZFR 


380  MAT  Qr«ZER 
390  MAT  Oi«ZER 

-AQ_Q _ MAT  Fr-mN 

410  MAT  Si*ZER 
420  MAT  Hi*ZER 
430  Nn«13Q«PI 


440  MAT  2r«(N0) 

450  MAT  Zi-ZER 

46,0 _ MAT  Cr«  7FP 


470  MAT  Gi*ZER  - 

480  Nrt*50 

■A9D _ PRINT  PAGF _  _ 

500  INPUT  -ENTER  NUMBER  OF  ROOTS  TO  BE  FOUND  (50)-, Nrt 
Xd  d* . 2 

■En«Qqr«l _ _ _ _ 

Qq i»Fc"F***0 


51  0 
Z2Ji 
530 
540 


E0-8 .854E-12 

-550 - R  g«6378Q_QS»4/3 _ 

560  U0-PI#4E-7 

570  C»1/SQR(E0*U0) 

_580_  INP_UT_-ENXERJEREQUENCY  _<MHzJ_  -  ,  F 
590  U*PI»F*2E6 
600  Lfs«C/(F*lE6> 

61Q Kf*-2«PI/Lfs 


620  FOR  L«1  TO  10 

630  DISP- “ENTER  DIELECTRIC  CONSTANT  OF  LAYER 
4 0 _ I N P UT  Er < L) 

650  IF  Er ( L ) <  >-l  THEN  680 
660  IF  L-l  THEN  630 
670  GOTO  730 

680  DISP  “ENTER  CONDUCTIVITY  OF  LAYER  “;L; 
690  INPUT  Si <  L )  9 

-TOfl_»D I "ENTER— TH LCKNES^nr  t  *vrP  .  . ,  . 


* ;L) “  (-1  TO  STOP)1 
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710  INPUT  Hi  <L) 

720  NEXT  L 

7 P14»WYM _ 

740  INPUT  “DO  YOU  WANT  HARD  COPY?  <Y/N)H,P1* 

750  IF  Pl*-'^"  THEN  PRINTER  IS  7 , 6 , WIDTH ( 1 32 ) 

.760  NLa-L^l  ...  _  —  _  _ _ 

770  PRINT  "LAYER  *  M  j  0  j  *  ;  Eps-"  j  1  >TAB  (30  )  j  "CONDUCTIVITY-- j  0 
780  FOR  L-l  TO  Nla-1 

_75LQ _ ERXbULJLLAY-ER-t  M;1;M  :  Epwl^fic (U>  jJEABX211  1 ;  wCONDiUCXXVI-TY>"  ^BiiL)  ;_TAB<  601^J1IHICKNESS«JLlHUUx^*»T*p*  h 

800  NEXT  L 

810  PRINT  "LAYER  *  "jNIjj":  Ep*»"  ;  Er  <Nla  >  j  TABOO  )  ;  “CONDUCTIVITY*"  >Si (Nla ) 

B2.Q _ Al-CKf**R«ZZl^Liy31 _ _ 

830  DISP  "  WAIT  -  SURFACE  IMPEDANCE  CALCULATIONS" 

840  FOR  L*1  TO  Nla 

-BM _ Cr(L)— Er(L)KE0ttU*WiOUI _ 

360  Gi<L)*Si(L)*W#U0 

870  PRINT  "Gartfn  " j L , Cr (L ) ,Gi < L ) 

880  NEXT  L  _  _ _ _ _ 

890  Niter— 1 

900  FOR  N*1  TO  Nrt+1 

9.1.0  -JCALL  Zi  u  r  f  (  T*.r  (  N )  ,Tti(N)  ,Zr  ,Zi) _ _ 

920  Der *Zr /NO 
930  Dei*Zi/N0 

940  Qr  <  N ) *A1  *Dei  _  ..  _ _ _ 

950  Qi ( N ) *-Al *Der 

960  IF  Qr  < 1  )  >0  THEN  1000 

-2ZQ _ Ea&*1 _ _ _ 

980  Xdd-1000 

990  GOTO  1010 

U-Q  a  Q.  _N  EXT _ N _ _ _ 

1010  Ni ter*Ni ter +1 
1  020  Qr  *Qr  <  1 ) 

1 03  0  Qi «ai  ( 1  1 _ 

1040  Qnag*SQR(Qr*Qr+Qi#Qi) 

1050  Argq*l 80/PI*ATN ( Qi/Qr ) 

__1  060_Fn«l. _ _ _ _  _  _ _ 

1  070  IF  Q«ag  < *1  THEN  1100 
1080  Fn*2 

tlOO  PR  INT^^F*"*;  F  ;  "" MH^^Qr^^Qr  ^  "Qi*'*  ;Qi 

1110  PRINT  CHR*< 124 )4“Q"4CHR$( 124)6 ;  SQR < Qr *Qr+Qi*Qi ) >TAB< 30) J " <Q*“ j Argq , LIN< 1 ) 

1  1  7n  PRINT  CHR$  (12A14J1D  “ 4CHR*.tl 24 )  4  M*  liSQR.C  D  e  r  *  D  e  r +JD  •  i  «  D •  i  ) jJ.AB <30_)j  "<D«"  ;  1 80/P  I»ATN  ( Dei/Der  ?  _ 

1130  PRINT  LIN< 1 >; "NUMBER  OF  ITERATIONS  *  " jNiter jLINCl ) 

1140  IF  Fss*l  THEN  1450 

JL1_5.Q._  J  a  d  d  *F  c  *  0 _ _ _ _ _ _ 

1160  IF  Ar qq  < -30  THEN  1300 
1170  IF  Argq  >-19,2  THEN  1200 

__UL8Q _ Fn*l _ _ _ _ _ _ _ _ _ 

1190  GOTO  1300 

1200  IF  Qnag<*l  THEN  1300 

121Q— LFL  Qwag.)3_  THEN  -1.270 _ _ _ _ _ 

1220  IF  <  Q*ag >2.5)  AND  <Argq>-.01)  THEN  1270 

1230  CALL  Kutta<0,0>200>QrlQi/Yilr<l),Yili(l),T*r(l),T*i<l),i) 

1240  HALL  RtQ  <  Qr  ,Qi  ,  Atr  ,  Ati) _  _ _ _ _ _ _ _ _ 

1250  PRINT  "Tsr <g)-“ ;Atr , "Tsi (g )«“ > Ati 
1260  GOTO  1280 

1270  CALL  R 1 0  (Qr  ,Qi , Tsr  < 1 )  »Tsi ( 1  )  ) _ _ _ _ _ _ _ 

1280  PRINT  USING  1390  ;  "Tsr  <  " ,  1  Tsr  ( 1  )  ,  "Tsi  <  “  ,  1  ,Tsi  (1  ),  "Qr*  ",  Qr  <  1)  ,  “Qi*"  ,  Qi  ( 1  ) 

1290  Jadd* 1 

_1 300-  FOR  J*1  TO  Nr_t_ _ 


1310  Js*J+Jadd 

1320  DISP  "  ITERATION  *“;Niter;"  ROOT  ♦  ";Js 

1330  IF  Fn=2  THEN  1360 

1340  CALL  Kutta (0,  0 ,50 ,Qr  < Js) ,Qi< Js) ,Yilr < J) ,Yili(J) ,Tsr <Js) ,Tsi< Js> ,1 > 
1350  GOTO  1380 


qr ,Qqi ) 


1370  CALL  Kutta(0,0,50,Qqr ,Qqi,Yi2r<J) ,Yi2i<J> ,Tsr<Js) ,Tsi<Js> ,2) 
1380  PRINT  USING  1 390 j "Tsr C, Js ,Tsr < Js >, "Tsi (“, Js Tsi ( J* > 

- -  1*  DDDD  >  _ 


1400  NEXT  J 
1410  BEEP 


“Qr* 


Qr ( Js) 


JL420  QS*"_je _ _ _ _ 

1430  INPUT  "Iterate  Roots  Again  (Y/N)?*,Q* 
1440  IF  QS* "Y "  THEN  900 


N-umber.  Cin,  ktt) 


1460  Dd*Din*1000 


,Din  ,Dm»l  t_,_Nd_i. 


1470  Ns*Nrt+Jadd 

T4f?n  uie  2.0  *  i  CTd.fc  /  (?»PT  )  ) _ 

1490  W«*l 

1500  IF  Qnag ) 1  THEN  1520 
1510  CA1  1  ArtQr,Qi) 

1520  FOR  Kin*l  TO  Ndi 
1530  X*Al»Dd/R* 

-L5AO — IF  _X->  Xdd— THEN-  L66Q _ _ _ _ — 

1550  IF  Qnag  < 1  THEN  1620 

1560  CALL  Scurv(X,Qr>Qi,Wr ;Ui,Rns>Nfp ,her ) 
1 570— IF— Nf  o*2— THEN  160  0 -  - 


"Qi* 


Qi ( Js) 
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1580  N**"S«all  Cur va tur»-l " 

1590  GOTO  1680 

-16-00-  4-1-1— Cur v* t-ur-«^3-g - 

1610  GOTO  1680 

1620  CALL  P*«ri«*<X,Qr ,Qi ,Ur ,Wi ) 

■  1630  NU"Pnw »r  Stih" _ 

1640  Rn««10 
1650  GOTO  1680 


1670  N*-"Re*idu*  S*ri**" 
1680  CALL  Cab*(Wr ,Ui ,Un) 


1700  PRINT  USING  171 0 jDd/1 000 ,Wpl#N$, "♦TERMS*  " ,Rn« , "Err or*" 
1710  IMAGE  4X,DDDD.DDD>9X,SDDD.D>10X,18A>5X,8A,DD,4X>6A,D.DE 

-1720  Dd«Dd*D«ul.t _ _ _ 

1730  NEXT  Kin 
1740  Q**"N" 


Q*- 


1760  IF  Q**"Y"  THEN  1450 
1770  Q*«"N" 

-L2ao _INPJJT_Ht«LAi*__R 0 o±*_Again_  <Y/N$7»,Q* 
1790  IF  Q*»"Y"  THEN  900 
1800  END  ! 


,M*r 


1010  DEF  FNTanh  <X) 


1830  IF  ABS ( X ) >10  THEN  1900 
1840  Ex-EXP<X) 


-1850  Tanh- <Ex-l/Ex)/< EjidOZE  xA 
1860  RETURN  Tanh 
1870  D2x-EXP <2*X>-1 

_L880._Tanh-D2x/C2+D2x> _ 

1890  RETURN  Tanh 
1900  Tanh*l 

— L91.Q  .RETURN— Tanjbi _ 

1920  FNEND  ! 


-193_0_  SUB  .Cd i  g i d  <  A 1  ,  B1  ,  A2  ,B2  ,R  ,  I ) 
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1940  OA2*A2+B2*B2 
1950  IF  COO  THEN  2000 

1960  PRINT  LINO),  "ERROR  IN  SUBPROGRAM  Cdivid  .  " 
1970  PRINT  "DIVISOR  IS  ZERO . " , LIN <2 ) 

1900  PAUSE 

-1990  SUBEXXT _ _  .  .  _ 

20  0  0  R=< A2*A1+B2*B1  )/C 
2010  1=  ( A2-X-B1 -B2XA1  )/C 

-ZQ2Q.  SUBEND — !. - 

2030  SUB  Cexp ( A , B , R ,1) 

ZQ 40  R=EXP(A)*COS(B5 
2050  I=EXP  (A)x-SIN(B) 

2060  SUBEND  ! 


2070  SUB  Cabs(X,Y,Cabs) 

2080  X1=ABS<X) 

2090  Y 1  =ABS  (  Y  ) _ _ _ _ _ _ 

2100  IF  Xl<  >0  THEN  2130 
2110  Cabs=Yl 

.2120  SUBEXIT _ 

2130  IF  Y 1 <  >0  THEN  2160 
2140  Cabs=Xl 

2150  SUBEXIT _ 

2160  IF  X 1  > Y 1  THEN  Cab s=Xl  *SQR  ( 1  +<  Y1  /XI  )  A2 ) 
2170  IF  XI  <  =Y1  THEN  Cabs  =  Yl  x-SQR  (1  +  (X1/Y1)A2) 
21£. Q SUBEND J 

2190  SUB  Csqr  t(A,B,R,X) 

2200  IF_ (AOQ)_OR.  TJBOfll  THEN  .223.(1  „  _  _ 

2210  R  =  I  =  0 
2220  SUBEXIT 

.2,230  CALL  Cabs  (  A  ,.B  ,.Caba) _ 

2240  R=SQR< (ABS<A)+Cabs)x .5) 

2250  IF  A <  0  THEN  2280 

2260  I  =  B/(R+R) _ _ _ _ 

2270  SUBEXIT 

2280  IF  B <  0  THEN  I=-R 

2290  TF  FO=n  THEN  I=R _ 

2300  R=B/(I+I) 

2310  SUBEND  ! 


2320  SUB  Zsurf (Tsr  ,Tsi > Zr , Zi ) 

2330  COM  W,EQ,U0 ,Er<10) ,Si<10> ,Hi<10) ,Gr<10) ,Gi<10> ,Kfs,Re,Nrt,Nla 

234  0  COM  Tsr ( 1 01 ) ,Tsi ( 101 )  ,Acr < 10) , Aci <10 ) ,Qr  <101 ) ,Qi<10i >  ,Zr <10  ).i2iU0  ) 

2350  Al=<Kfs*Re/2)A<l/3)/Re 

2360  Lsr  =t<f  s+Al  x-Tsr 

■  2370  L  s  i  =  A 1  X-Te^j _  .  _ _ 

2380  Lsr 2=Lsr*Lsr-Lsi*Lsi 
2390  Lsi2-2*Lsr*Lsi 

■24DQ  E_DR  M=Nla  TQ  2  STEP  -1  _ _ _ _ _ _ _ _ 

2410  CALL  Csqrt <  Gr  <M-1 >+Lsr2,Gi <M-1 )+Lsi2,l)ir ,Uii > 

2420  CALL  Csqr t < Gr ( M ) +Lsr2 , Gi ( M ) +Lsi2 , Uwr , Uwi ) 

_243Q_CAL2  CdivULCUir  ^Ui^L^SiXM-U^ELOx-Er-CM-n x-W _  _  _ 

2440  IF  M=Nla  THEN  2480 
2450  Kwr-Zr  <  M) 

_ 24 6  0 _ K.w  i_=Z.i_<  MJ _ _ _ 

2470  GOTO  2490 

2480  CALL  Cdi vicK Uwr  ,  Uwi  , Si  < M )  , EG*Er  < M > xW , Kwr  , K wi  ) 

— 24.9JQ_.  A  t.= ENT  a  n  h  <  ULcx  Ki  <  Mr  1 )  )__ _  _ _ 

2500  Ai=T  AN ( Ui  ix*Hi  ( M-l  )  ) 

2510  CALL  Cdivid (Ar ,Ai ,1 , Ar*Ai ,Tr ,Ti ) 

3520  N  r =K  ir*TrH<  ii*Ti  +  Kwr _ _ 

2530  Ni=K irxTi+K iixTr+K wi 
2540  Dr=Kwr xTr-K wixTi+K ir 
_2550fcEi=KwrxTi+Kwixrr+Kii 
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2560  CALL  Cd i v i d ( Nr , Ni , Dr , Di , Ar , Ai ) 
2570  Zr<M-l >=Kir*Ar-Kii*Ai 
2580  Zi <  M-l  )  =K ir *Ai+K ii*Ar 
2590  NEXT  H 
2600  Zr=Zr ( 1 ) 

261 0  Zi=7i H  ) _ 

2620  SUBEND  ! 


^KerjL 


2640  COM  U.EO ,U0,Er(10) ,Si<10> ,Hi(10> ,Gr (10) ,Gi(10> , Kf s , Re , Nrt  > Nla 
2650  COM  Tsr  < 1 01 ) ,Tsi ( 1 01 ),Acr(10) , Aci < 1 0 > ,Qr < 1 01 > , Qi ( 1 01 > ,Zr ( 1 0 ) ,Zi < 1 0) 

-266 0-  Ec * F-i - S r - S i i * FLs*.  Q_ _ _ _ _  _  _ _ 

2670  Er=l E-6 

2680  CALL  Cd i v id < P I *X , 0 , 0 , 1 , Ar , Ai ) 

2690  r.ALL.Osqp.t  (  Ar  ,  Ai  ) _ _ _ 

2700  FOR  S*1  TO  Ns 


271  0  Qr2=Qr  <S)*Qr  (S)-Qi  <S)*Qi  <S> 

■272.Q,  Qi2=Zx-gr  C£l»Q.i  (3) _ 

2730  CALL  Ce xp  ( Tsi  ( S )  *X  , -Tsr  <  S )  *X ,  Ar  ,  Ai  ) 

2740  CALL  Cd i v id ( Ar , Ai , Tsr < S ) -Qr 2 , Tsi ( S ) -Qi2 , Cr , Ci ) 
2.7-5Q-F.r=Fr+Cr _ 

2760  Fi=Fi+Ci 

2770  IF  S<10  THEN  2830 

27J3il .  Err =ABSlC  (Fr -Sr  )./Sr  ) _  _ _ _ _ _ 

2790  Eri=ABS( ( Fi-Si i ) /Si i  ) 


2800  IF  Err  >Er  THEN  2830 
-2S..1-0— IF...  Er  i  >Er  -T.HEN-2S3Q 

2820  Fs=l 
2830  Sr=Fr 


J2BASL  SH=F-i _ _ 

2850  IF  Fs*1  THEN  2870 
2860  NEXT  S 

■287.0-Hen=MAX  ( Err  ,  Er  i  )_ 

2880  Wr=Br*Sr-Bi*Sii 
2890  Wi=Br*Si  i+Bi-xSr 

29QQ-RnssrS-l _ 

291 0  SUBEND  ' 


.25.211  SUB  .  1<  u XX a  X  Ar  ,  Ai  ,  N  ,  hr  ,  Hi  ^alr^..Yn.ti...,  Y r -,Y.i.,Fn> 


2930  RAD 
2940  DIM  Kr ( 3 ) ,Ki(3) 
255.0  Hr  =XIlr  -  An.)  /  N _ 


2960  Hi“< Bi-Ai ) /N 
2970  Eps=l E-6 


2990  Xi=Ai 
3000  Hhr=Hr/2 

_30_LQ _Hhi=Hi/2 _ 

3020  Ysvr=Yr=Ynfr 
3030  Ysv i=Y i=Yn  t i 
313  40  Xsvr  -Xr 


3050  Xswi=Xi 
3060  FOR  L=0  TO  3 


:i,Fr,Fi) 


3080  Kr(L)-Hr*Fr-Hi*Fi 
3090  Ki<L>=Hr*Fi+Hi*Fr 

31  Q  Q  ON  L±JLXQXa  3110,3110,3160  ^.20JL 
3110  Xr=Xsvr+Hhr 
3120  Xi=Xs v i+Hh i 


-313 CL_ Y_sv.r=Y r+Kr (L)» .5 . 
3140  Ysui=Yi +Ki(L)*.5 
3150  GOTO  3200 

3160  Xr=Xs v r±Hr _ 

3170  Xi=Xsvi+Hi 
3180  Ysvr=Yr+Kr <L> 

.3  \2  °___Ys y  i  =  Y  i  +K  i  <  L ) 
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320  0  NEXT  L 

3210  Yr=Yr+(Kr  <  0  )  +2*  ( Kr  ( 1  >+Kr(2)  >*H<r<3>  )/6 
3220  Y  i = Y  i  +  ( K  i (0) +2* ( K i ( 1 ) +K i C 2 ) )+Ki(3) >/6 
3230  IF  Xr>Br-Eps  THEN  SUBEXIT 
3240  Ysur=Yr 


-3250., 

3260 

3270 

Y  i  - 

GOTO  3040 

SUBEND  ! 

3280  SUB  F  u  n  c  <  F  n  >Ysvr>Ysvi,Xr>Xi>Fr>Fi) 

329  0  Ar=Xr*Xr-Xi*Xi 

330  ft  A*i  =2*Xr*Xi _  _  .  .. 

3310 

3320 

33 3  fl 

IF  Fn=2  THEN  3340 

CALL  Cdi y  id  ( 1 ,0  > Ysyr-Ar ,Ysvi-Ai,Fr  > F i ) 

81 IRFXTT 

334  0  Br=Ysvr*Ar-Ysv i*Ai 

3350  Bi=Ysvr*Ai+Ysvi*Ar 

3360  CALL  Cdivid < 1  .  0  .  1-Br  .-Bi .Fr  .Fi > 

3370 

SUBEND  ! 

3380 

SUB  npnupr(N,Xr;XiJYr;Yi) 

3390 

3400 

3410 

IF  NOD  THEN  3430 

Yr=l 

Yi  =  0 

3420 

3430 

3440 

SUBEXIT 

Yr=Ar=Xr 

YisrAisXi 

3450 

3460 

3470 

IF  N=-l  THEN  3570 

IF  N*1  THEN  SUBEXI1 

FnR  1  =1  TO  ABS(N>-1 

3480 

3490 

3500 

Zr=Ar*Xr-Ai*Xi 

Zi=Ar  -x-Xi+Ai  -x-Xr 

Ar  —7  r 

3510 

3520 

3530 

Ai=Zi 

NEXT  L 

TF  N<  n  THEN  3570   . 

3540 

3550 

3560 

Yr  =Ar 

Yi=Ai 

SUBEXIT  . . 

3570 

3530 

CALL  Cdivid <1 , 0 ,Ar ,Ai , Yr , Yi> 

SUBEND  ! 

3590 
3600 
361  0 

SUB  Rt0(X,Y,T0r ,TQi> 

Ar=Ai“Br=Bi=Cr=Ci=0 

CALL  Cp  ower  <2.^X  .Y  .Ar  ^  AL> 

3620 

3630 

3640 

CALL  Cpower(-i ,X,Y,Br,Bi) 

Ar=Ar +  ,  5*Br 

Ai=Ai+  5-x-Bi 

3650 

3660 

3670 

CALL  Cpower<-4,X,Y,Br ,Bi) 

Ar=Ar+Br/8 

A*i  =A  i  4-B  i  /R 

3680 

3690 

37nn 

CALL  C’power(-7,X>Y,Br,Bi) 

Ar=Ar+5*Br/32 

A  i  =  A  i  +  5  x  B  i/32 

371  0 
3720 

CALL  Cpower (-1 0 ,X,Y,Br ,Bi> 

Ar=Ar  +  l  l  -xBr/32 

AiaAi-f  1  1  X-B i  /37 _ 

3740 

3750 

3760 

CALL  CpowerC-lZ.XjYjBr ,Bi) 

Ar=Ar+539*Br/512 

Ai=Ai  +  539*x- Bi/512 

3770 

3780 

3790 

CALL  Cpower (3,X,Y>Cr,Ci) 

Cr=-l-Cr*4/3 

Ci=-Ci *4/3 

3800 
381  0 
3820 

CALL  Cp ower < -3 , X , Y , Br , Bi ) 

Cr=Cr-7*Br/12 

Ci=Cj.x7*&L/l2 
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3830  CALL  Cpower ( -6 , X , Y , Br , Bi ) 

3840  Cr=Cr-31 *Br/48 
3850  Ci*Ci-31*Bi/48 
3860  CALL  Cpower  <-9,X,Y,Br,Bi> 

3870  Cr=Cr-397*Br/288 

388(1  Ci  =Ci. -397xBi/2B  8 _ 

3890  CALL  Cexp ( Cr , Ci , Br ,Bi ) 

3900  CALL  Cpower (2, X,Y,Cr ,Ci> 

3910  TO e_=Ae_+  2*B c_* Ci  +9x- B i_x.Cc _ 

3920  T0i=Ai+2*Bi*Ci-2*Br*Cr 
3930  SUBEND  ! 

3940  SUB  Azer o <Z ,  R  > 

3950  R=l+5/(48xZA2>-5/(36xZA4)+77125/(82944*ZA6>-l  08056875/ ( 6967296x-ZA8  > 

-3.960 - Rg-XA-C2/3 X xXR-H  69375596875/  l33.443Q.2Q8xZ21.0_) J _ _ 

3970  SUBEND  ! 

-39 a 0-SUB  Az.ec  op  (Z^fU _  _ _ 

3990  R=l-7/(48xZA2)+35/(288xZA4)-181223/(207360xZA6) 

4000  R=ZA(2/3)*(R+l 8683371 /(1244160xZAS)-911 45884361 / (191102976xZAl0) > 
4010  S1IRFNP  i 


4020  SUB  Pser ies (X , Qr , Qi ,Wr ,Wi ) 

403.Q  .CO«_U>Ea.,U0  *£eX10.)  ,Sitl.0)  ,tliC10>  ^GctlQ)  ,Gi(10)  ^Kf  s  ,  Re  ,  Nr  T  ,  Nla 
4  040  CDH  Tsr ( lOi > ,Tsi (1 01 ) ,Acr (10) , Aci < 1 0 > ,Qr < 1 0 1 > , Qi < 1 01 > ,Zr ( 1 0 > ,Zi ( 1 0 > 
4050  Wr =Wi=0 

-4Q60  CAM  r.P»p(0JPT/4)Ar,A-i'> _ _ 

4070  Br=<Qr*Ar-Q-*Ai)*SQR (X) 

4080  Bi= ( Qr xAi+Qix*" ) xSQR ( X) 

—40-9  0 — EQR_  M-O-  XQ—LO _ _ _ 

4100  CALL  Cp ower (M ,Br ,Bi , Ar , Ai ) 

4110  Dr=Acr  (M)  x-Ar-Aci  (M)*Ai 


4130  Ur=Wr+Dr 
4140  Ui=Ui+Di 

-A  150-NEXX  -M _ _ 

4160  SUBEND  ! 

4170  SUB  ArOOr.ni  1 _ 

4130  COM  W , E0 ,U0 , Er ( 1 0 > , Si <1 0 > , Hi ( 1 0 ) ,Gr ( 1 0 ) , Gi ( 1 0 > , Kf s , Re , Nr T , Nla 
4190  COM  Tsr (101 ) ,Tsi (101 ) , Acr ( 1 0 ) , Aci ( 1 0 ) , Qr ( 1 01 > ,Qi < 1 01 > ,Zr ( 1 0 > ,Zi < 1 0 > 

-420  0  P 2=&QR_CE U _ 

4210  Acr ( 0 1 =1 
4220  Aci(0>=0 

4230  Arr  (1  1  =  0 _ _ _ _  _  _ 

4240  Aci ( 1 >=-P2 
4250  Acr(2>=-2 

4260  Ar -i  (91  =  0 _ 

4270  CALL  Cp ower ( -3 , Qr , Qi , A3r , A3i ) 

4280  CALL  Cp ower ( -6 , Qr , Qi , A6r , A6i ) 

—4290  CAU — Cpn*ec.l=9.,rjr- ,  Qi  ^A9c-^A9iJ. _ _ _ 

4300  Acr (3>=-(A3i/4)*P2 
4310  Aci (3>  =  ( 1+A3r/4)*P2 

4390  Arr  (41=4/3 -X-(1+A3r /2) _ _ _ 

4330  Aci  ( 4>=4/3*A3i/2 
4340  Acr(5)=P2/4*(3xA3i/4) 

— 43.50 _ A  Ci  (.5 )  =- ( P2/4 )  *  ( 1  +3  *A3r/4 > : _ 

4360  Acr ( 6 ) =- ( 8/ 1 5 ) * < 1 +A3r +7*A6r /32 ) 

4370  Aci(6)=-(8/15)x(A3i+7*A6i/32) 

—4330  Acr_(_? )  =— (P5/6)x( 5*A3i /4+27xA6i  /32.) _ _ _ _ 

4390  Aci  (7)=P2/6*(  l+5x-A3r/4+27xA6r/32) 

440  0  Acr ( 8  >  =  1 6/1 05 x ( 1  +3*  A3r /2+27* A6r  /32  ) 

441  0  ftri  (81  =  16/1  05*(3xA3i  /9+97*66i  /39  1 _ _ _ _ _ 

4420  Acr  (  9)=P2/24*(7*A3i/4+5*A6i/4+21  *A9i/64  ) 

4430  Aci (9)=-(P2/24)x(l+7xA3r/4+5xA6r/4+21xA9r/64) 

,  4440  Acr  (  1  0)=-  (32/945 +64-xA3r/945+l  1  *A6r /I  89+7*A9r/270  ) 
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445 Q  Aci  < 10)=-( 64* A3 i/945+11 *A6i/ 1 89+7* A9 i/270 > 

4460  SUBEND  ! 

4470  SUB  Scurv(X,Qr ,  Qi  ,  Wr  ,  I4i  ,  Sn s  ,  Nfp  ,  Her  > 

4430  Wr =Wi=Q 

-4490 — Ac  a  dr  *dr_=£UJfc&i - 

4500  Ai=2*Qr*Qi 
4510  Pr=-X*Ai 

-4520...  P.i.=XxAr - — 

4530  PMag=SQR (Pr*Pr+Pi*Pi ) 

4540  Psr=Prx*Pr-Pi*Pi 

-4  5  5  Q  P  s  i  = -2.x-  ELc-*.Ri - — -  — 

4560  IF  <  Pwag  >12)  OR  <Qr<0)  THEN  4600 
4570  CALL  Fp 1 < P r  , P i , Fr , Fi , Sns , Her > 

45R0  NFp  =  1 _ _ 

4590  GOTO  4620 

4600  CALL  Fp2<Pr,Pi>Fr ^Fi.Sns.Mer) 

_A61Q,  Nfp^2 - - - 

4620  CALL  Csqr  t  C  Pr *P  I  > P  i *P  I ,  P2r  , P2i  ) 

4630  CALL  Cp o wer ( -3 , Qr > Qi , A3r > A3i ) 

4640  CALL  Cp  oner  ( -6  „ Qr., QL,  A6r  »  A6i.) -  - — 

4650  Br«l+P2i-Fr-2*Pr*Fr+2*Pi*F i 
4660  Bi«-<P2r+Fi+2*Pr*Fi+2*Pi*Fr ) 

_Ah2SLMnsFr+  .25*  ( A3r*Br.rA3i*Bi  ) . .  . . . . 

4600  Ui=Fi+,25*<A3r*Bi+A3i*Br > 

469  0  Br=l-P2r*Pi+P2i~P2i*Pr-2*Pr+5/6*Psr+.5*Psr*Fr-Fr- ,  5*Psi*Fi 
4700  B  l  ~  P  2  r  *  P  r-  -P-2c  -  P  2  -i  *  P  A  -  2  x  P  i  +  5/  6  *  P  jsi.  ■ +  5  *  P  sr  .*  F_i -F  i+.  ^5  *  P  s  i  *Fr 
4710  Ur=Wr+ , 25*  <  A6r #Br-A6i*Bi ) 

4720  Wi=Ui+ . 25* (  A6r*Bi+A6i*Br  ) 

AZ3Q_  .SUBEND.  1 _ _ _ _  _ _ _ _  - 

4740  SUB  Fpl (Pr ,Pi,Fr ,Fi,Sns, Her ) 

.4750  Cr=Ci-Fr=Fi=Fs-Q _ _ _ _ - _ _ 

4760  Er  =  l E-4 

4770  CALL  Csqr t < Pr , P i , Ar , A i ) 

_A78Q  _Tr=-Pr _ _ ... _ _ _  _  _ 

4790  Ti—Pi 
4800  Gr=Br=-Ai 

. 48.1.0-  Gi-Bi^Ar _ _ 

4320  FOR  N*1  TO  40 

4830  F 1 =- < 1 /N  >  *  < <2*N~1 >/<2*N+l ) ) 

.Aa4.Q.  Cr=Fl*(Br*Tr.-Bi.*Ti)  _  _  _ _ _ _  „ 

4350  Ci=Fl *<Br*Ti+Bi*Tr ) 

4860  Br =Cr 

4870  Bi=Ci _ _ _ ... 

4880  Gr=Gr+Br 
4890  Gi=Gi+Bi 

_490  0  IF  N<3  THEN  4960  .  „  _ .. _ _ _ 

4910  Err=ABS( (Fr-Gr)/Fr) 

4920  Eri=ABSC <Fi«Gi)/Fi) 

4930  IF  Err) Er_THEfl_496Q _ _ _ _  _ 

4940  IF  Eri>Er  THEN  4960 
4950  Fs=l 
4960.Fr^Gr 
4970  F i =G i 

4730  IF  F*«l  THEN  5000 

4990  NEXT  N _ _ 

5000  Mer=MAX <Err ;Eri ) 

5010  Sns=N-l 

—5.0  20  Cr  =  l.r 2/SDR  (PI  )»Fr _ - _ 

5030  Ci»-<2/SQR  <  PI ) ) *Fi 
5040  CALL  Cexp < -Pr i , Br , Bi ) 

5050  Gr^Erx-Cr-Bix-Ci _ 

5060  Gi=Br*Ci+Bi*Cr 

5070  CALL  Csqrt <Pr*PI ,Pi*P I ,Ar , Ai) 

_  500  0  Fr^j+ArxGi+AixGr _ _ _  _  _ 
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5090  Fi=Ai*Gi-Ar*Gr 
5100  SUBEND  ! 


5110  SUB  Fp2(Pr,Pi>Fr>Fi>Sns ,Mer ) 
5120  Gr=Gi=Fr~Fi=Fs=s0 
-5130 . Fr  =  l E-4 _ 


5140  IF  ATNCPi/Pr  >  <  0  THEN  5190 
5150  CALL  Csqrt <PI*Pr ,PI*Pi ,Cr ,Ci> 
U-6.fl-.CAiJ  Cp*p(-Pr  >  


5170  Gr-2*Dr*Ci+2*Di*Cr 
5180  Gi=2*Di*Ci-2*Dr*Cr 
- .5JJ9 0_CAU —  Cri-i-v irLCL^Q. ,  -  2  •*  P n_,  -2*P  i ,  Ar ■.  *  A i  )  „ 
5200  Tr=A r 
5210  Ti=Ai 
522Q. Gr-Gr-fTr 


5230 

5240 

5260 

5270 

5280 


Gi=Gi+Ti 
FOR  M=2  TO 


Bi--  <  2*M-1  >*< Ar*Ti+Ai*Tr ) 
Tr=Br 
XiJsBi 


5290 
530  0 


Gr~Gr+Tr 
Gi=Gi+T i 

XHEN- 


3Z0- 


5320  Err=ABS((Fi — Gr)/Fr) 
5330  Er i-ABS ( (Fi-Gi ) /Fi ) 


- 

5350 

-xt — t.ncvfc.r — u-ttUNf— zo/.u - - - - 

IF  ErDEr  THEN  5370 

5360 

Fs=l 

537  0 

.  F.r.=Gr _ 

5380 

Fi=Gi 

539  0 

IF  Fs-1  THEN  5410 

-5AQ-Q- 

-NE-X-T-  H - 

5410 

Ner=MAX (Err , Er i  ) 

5420 

Sns-N-1 

J5-4-3A 

-SUBEND . 
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Appendix  B 

USEFUL  APPROXIMATION  FOR  THE  GROUND-WAVE  ATTENUATION 


1 


Solutions  for  the  ground-wave  attenuation  function  can  be  found  by  using  one  of  three  tech¬ 
niques.  Two  of  these  are  approximations  valid  at  short  numerical  distances  while  the  third  uses  a 
residue  series  summation  that  is  technically  accurate  at  any  range,  but  which  becomes  increasingly 
inefficient  (in  the  number  of  terms  required)  as  the  distance  between  source  and  receiver  is  decreased. 

Calculation  of  the  attenuation  function  using  the  residue  series  representation  requires  knowledge 
of  the  poles  ts,  which  are  roots  of  the  equation: 

^P--qw(t)=0,  (Bl) 

at 

where 

w(t)  =  \TK{Bi{t)  —  iAi(t) , 
in  which  Ai(t)  and  Bi(t)  are  Airy  functions  [Bl]. 


An  efficient  procedure  for  solution  of  Eq.  (Bl)  has  been  described  by  Hill  and  Wait  [B2].  A 
summary  of  this  procedure  and  the  regions  of  validity  are  presented  in  Fig.  Bl.  In  the  first  region, 
denoted  by  “1,”  the  differential  equation 


dts_  _  l 
dq  ~  ts  —  q2  ’ 


(B2) 


is  integrated  numerically  using  a  complex,  Runge-Kutta  technique  (Appendix  A  provides  a  listing  of 
the  numerical  techniques  employed  to  solve  for  the  attenuation  function)  starting  from  the  initial  con- 
ditions: 


=  o  =  ',(0)  =  0i'se 


l  7T 

3 


in  which  the  set  of  values  a 's  (s  =  1  to  oo)  are  the  real  zeroes  of  Ai'(-a). 


(B3) 


In  the  second  region  identified  by  “II”  in  Fig.  Bl,  the  transformation  Q  =  q  1  is  made  to 
obtain  the  differential  equation: 

dts  i 

—  =  - (B4) 

dQ  1  -Q2ts 

This  equation  is  integrated  using  the  same  numerical  technique  starting  from  the  initial  condi¬ 


tions: 


ls  I  Q 


=  0  =  **(°°)  =  ase 


l  7 T 

3 


(B5) 


in  which  the  values  as  (s  =  1  to  oo)  are  the  real  zeroes  of  Ai(— a). 


In  the  region  of  Fig.  Bl  in  which  the  ground  constants  produce  |  q  |  >  1  and  phase  of  q  >  —30°, 
an  additional  root  t0  represents  a  trapped  surface  wave  mode.  Depending  on  the  exact  value  of  q, 
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Phase  of  q 


II. 


dtr 


dQ  (l  -  Q\ ) 


t,  \q  =  0  = 


where  Q  —  q 


_  „-i 


in.  ,0  =  ,2  +  J_  +  _L  +  ^  +  _^  +  _mT  + 


IV. 


dtc 


2  q  8  q4  32q7  32q10  5l2q13 

1 


dq  (ts  -  q1) 


*0  I  q  =  0  —  rl(0) 


Figure  B1  —  Methods  used  for  solution  of  the  roots  ts  and  their  regions  of  validity 


two  methods  are  used  to  obtain  this  root.  The  first  (given  as  “IH”  in  Fig.  Bl)  uses  the  asymptotic 
expansion  [B3]: 


~  q2 

—  I2q2exp 


+  +  — h-  + 

2  q  8  q4 


4  q3 


-  1 


+ 


11 

32(? 10 
31 

48<76 


+ 


539 
512^13 
397 
288 q9 


(B6) 


The  second  method  (given  as  “IV”  in  Fig.  Bl)  uses  equation  (B2)  in  the  region  shown  for  the 
solution  of  the  root  t0  only,  with  the  initial  condition  t0\q  =0  =  *i(0).  All  other  roots  within  the 
applicable  region  are  found  using  Eq.  (B4). 


Two  approximations  are  useful  for  calculation  of  the  ground  wave  attenuation  function  at  short 
numerical  distances  x .  The  first  of  these  is  a  modification  to  the  flat  earth  attenuation  function  for 
spherical  earth  and  is  given  [B2]  by: 


W  =F(p)  +  ~~ 
Aq 


1  —  i  Vrp  -  (1  +  2p)F(p) 


+ 


1 

4  q6 


1  —  i  xfnp  (1  —  p  )  —  2p 


4- 


F(p) 


(B7) 
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where 

F(p )  =  1  —  i  yfirp  e('~p^ erfcii'fp) ,  (B8) 

in  which  p  =  ixq 2  and  erfc  is  the  complementary  error  function  [Bl], 

The  approximation  is  useful  at  short  numerical  distance  x<0.2  and  for  values  of  |  q  \  >  1  as 
shown  in  Fig.  B2.  Typically,  these  conditions  can  apply  to  a  uniform  lossy  surface  as  well  as  certain 
layered  cases. 


Nurrerical  Distance  x 


x  <0.2 


x  >0.2 


I.  Power  Series 

II.  Small  Curvature 

III.  Residue  Series 


Figure  B2  —  Methods  used  for  calculation  of  the  ground  wave 
and  the  regions  over  which  they  are  most  suitable 


The  second  short  distance  approximation  is  a  power  series  due  to  Bremmer  [B4]  given  by: 

j  7T 

w  =  E  An(e  4  W*)". 
n  =0 

in  which  the  first  five  coefficients  are: 

^o=l, 

A  j  =  —  i  V  it  , 

A2  =  -2, 

1 


(B9) 


A%  =  i  "V  7T 


Aa  = 


1  + 


1  + 


4  qi 
1 


and 


2  q3 


Additional  terms  through  A 10  are  available  [B2]  and,  while  they  are  not  included  here,  have  been  used 
in  the  calculations  presented  in  this  report.  The  power  series  approximation  is  most  useful  for  a 
numerical  distance  x  less  than  approximately  0.2  and  for  a  surface  characterized  by  |  q  \  <  1  such  as 
unlayered  seawater  or  thin  ice  over  seawater  at  low  frequencies. 
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